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C _~) • The Sturm— Liouville eigenvalue equation for eigenmodes of the radial oscillations is de- 

£N| ' termined for spherically symmetric perfect fluid configurations in spacetimes with a 

nonzero cosmological constant and applied in the cases of configurations with uniform 
distribution of energy density and polytropic spheres. It is shown that a repulsive cosmo- 
logical constant rises the critical adiabatic index and decreases the critical radius under 
which the dynamical instability occurs. 
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1. Introduction 

t-h ; 

. The internal Schwarzschild spacetimes with nonzero cosmological constant (A ^ 0) 

and uniform distribution of energy density were given for star-like configurations 1 
■ and extended to more general situations. The polytropic and adiabatic spheres were 

t^- , preliminary treated and compared, 3,4 neutron star models with regions of nuclear 

matter described by different relativistic equations of state that are matched were 
also treated. 5 Their stability can be considered on energetic grounds 6 or it can be 



treated in dynamical way. 7 Here we determine the dynamical stability conditions 
for the uniform density and polytropic spheres using the approach of Misner et al. 8 



2. Sturm Liouville Equation 



x. 

Interior of a spherically symmetric configuration is described (in standard 
Schwarzschild coordinates) by the line element 

d s 2 = -e 2 * dt 2 + e 2 * dr 2 + r 2 (d0 2 + sin 2 0d^ 2 ) , (1) 

with metric coefficients taken in the general form 

* = *(r,i), $ = $(r,i). (2) 

The perfect fluid distribution is given by energy density and pressure profiles p(r, t) 
and p(r,t). The static equilibrium configuration is given by $o( r )> ^ / o( r ) ) Po( r )) 
Po(r) satisfying the Euler equations. The radially pulsating configuration is then 
determined by 

q(r,t) = q Q (r) + 8q(r,t), Sq = (5$, <S* , Sp, 8p, Sn) , (3) 
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where n is the number density. The pulsation is given by the radial displacement 
£ = t). The Euler perturbations Sq are related to the Lagrangian perturbations 
measured by an observer who moves with the pulsating fluid by the relation 

Aq(r, t) = q{r + £(r, t), t) - q {r) « Sq + fa . (4) 

Introducing a new variable ( = r 2 e _ *°£, the radial pulsations are governed by 

W( = (PC')' + QC (5) 

with the functions W(r), -P(r), Q(r) determined for the equilibrium configuration 



\V= ( Po +Po)^e 3W \ ((>) 
P^ 1Po ±e*° + ™\ (7) 



r 

,*o+34>n 



Q = e 



(Po) 2 1 4p /8ttG ,\e 2 * 



- (po +Po) — - A 



(8) 



_Po+Por z r A \ c 

The linear stability analysis can be realized by the standard assumption of the 
displacement decomposition 

CM) = CMe iwt . (9) 

Then the dynamic equation reduces to the Sturm-Liouville equation and the related 
boundary conditions in the form 

{PC)' + {Q + " 2 w)t = o, (10) 

C e*° 

— finite as r — > , "fPo — ~ * as r~^R. (11) 

The Sturm-Liouville equation (|10[) and the boundary conditions determine eigen- 
frequencies uij and corresponding eigenmodes Ci( r )> where i — 1,2, ... ,n. The eigen- 
value Sturm-Liouville (SL) problem can be expressed in the variational form 8 as 
the extremal values of 

= (12) 

f Q R W( 2 dr 

determine the eigenfrequencies u>i and the corresponding functions Q(r) are the 
eigenfunctions that have to satisfy the orthogonality relation 



e 3 * 



°~*°(po + Po)r 2 ^ 4) C (, ' ) dr = 0. (13) 



3. Results and Conclusions 

We applied the Sturm-Liouville approach to spheres with uniform energy density, 1 
and polytropic spheres. 3 The case of uniform spheres can be properly taken as a test 
bed of the dynamical instability problem — although these solutions of the Einstein 
equations are of rather artificial character, they reflect quite well the basic properties 
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Figure 1. Dependence of the critical value of adiabatic index 7 C on sphere radius R. Full curve: 
vanishing cosmological constant A = 0; then 7 C diverges as R — * 9r g /8 from above. Dashed curve: 
positive cosmological constant A = 0.1, the point of divergence is shifted to 1.18421 > 9/8, and 
7c A>0 > 7cA=o- Dashed- dotted curve: negative cosmological constant A = —0.1, the point of 
divergence is shifted to 1.07143 < 9/8, and 7 c ,a<o < 7c,A=o- 



of very compact objects. 9 The dependence of the critical value of adiabatic index 
7 = (d\np/dlnn) s — (n/p)(Ap/ An) on configuration radius R is for the uniform 
case given in Fig. [T] (see also 12 ). The polytropic case is treated in details by Stuchlik 
and Hledik. 10 
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